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SYNOP SI S 

IN PLANE VIBPJiTION OP LAYERED RINGS ON 
PERIODIC RADIAL SUPPORTS USING 
FINITE ELH^IENT METHOD 

A Thesis Submitted 

In P artial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 
by 

J.V.SASTRY 
to the 

DEPARTMENT OF MECHANICAL ENGINEERING 
INDIAN INSTITUTE OP TECHNOLOGY, KANPUR 

JULY, 1936 

The thesis presents the analysis of the in-plane 
vibration of layered rings with periodic supports using 
finite element approach. 

MaMng use of wave propagation theory many 
researchers have found exact harmonic solutions to many 
problems. But the exact harmonic solution becomes very 
difficult to obtain, when anisotropy of the material is 
considered. In all such cases, thv. present P EM analysis 
would provide accurate propagation constants, natural 
frequencies and corresponding mode shapes. 

In the present workj^ a finite element cornputer 
programme in PORTRAN-IV is developed and is iir^lOTented 



(xiv) 


on DEC- 1090 system to get frequencies of a periodically 
supported beam and are found to compare well with the 
results presented by earlier x* 7 orKers. 

The finite element approach is also used to 
get frequencies and mode shapes of a two layered elastic 
ring with the following support conditions. 

Mono coupled 

i) The supports prevent both radial and tangential 
displacements (u and v respectively) resulting 
in a mono coupled system having au/fe©--.- as the 
only coupling coordinate between one bay and 
next bay of the ring« 

Si-coupled 

ii) The supports Prevent only u, the radial dis- 
placesment at the supports, thus mahing the 
system a Bi-coupled having v and su/^© as the 
coupling coordinates between one bay and the 
next bay of the ring and the natural frequencies 
are found to be within 0,5^ a ccuracny compared 

to those presented by Reddy [9] who used wave 
approach, 

Finally, the present method is applied to 
obtain freqaencies and mode shapes of a three layered 
orthotropic ring and these results are tabulated. 



CHAPTER - I 


INTRODUCTION 

1,1 Intjroducbion an<3 Literature Review? 

Many structures like Bridges/ Aircraft fuselage 
etc., contain a series of identical structural elements 
linked by identical junctions. These types of struc- 
tures, with a regularly repeating section, having a 
combination of support conditions are referred to as 
periodic structures. 

It is often required to study the dynamic 
behaviour of these structures for fluctuating pressure 
fields, frequently it is also necessary to assess 'their 
vibration inodes and response characteristics. Here it 
is interesting to note that, instead of finding -the 
normal modes and natural frequencies of -the structure, 
it is preferable to describe the behaviour of the periodic 
system by -the proper-ties of the characteristic waves, 
as this gives a better physical understanding of the 
syston behaviour. Each of these characteristic waves 
is. associated wi-th a propagation constant. Also> 
there is a constant ratio between -the amplitude of 
motion in one bay and that at the corresponding point 
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in the adjacent bay. The ratio of the ainplitudes is 
equal to e^, where n is the propagation constant, which 
may be in general a complex quantity. It's values 
always occur in positive and negative pairs correspond- 
ing to identical but opposite going waves. The real 
part of M is known as 'Decay Constant' and the imaginary 
part is called the 'Phase Constant', 

Wave propagation will depend on whether for the 
particular wave, the propagation constant is purely 
real, purely imaginaaq^ or complex. If the propagation 
constant is. purely real, the wave is non-propagating. 

In that case e^ represents ^n exponential decay or 
growth of displacement from bay to bay, with correspond- 
ing displacements of adjacent bays being in phase. If 
propagation constant is purely imaginary (say i^ ) then 
the factor e^ represents a simple phase difference of 
between the corresponding displacements of adjacent 
bays, but the amplitudes of the displacenents are the 
same. If propagation constant u has both non- zero real 
and imaginary parts, imaginary part being 0, e^ then 
represents exponential decay or growth of displacement 
but with corresponding displacements in adjacent bays 
being out of phase by an angle 0. 

Wave motion in simple periodic systems has 
been studied for nearly ^10 years, as Brillouin [l] 
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electrical engineers have developed the studies over 
the years in relation to Crystals, optics, electrical 
transmission lines, etc. W-ave motion in engineering 
periodic structures (consisting of beams, plates etc.) 
has been investigated in early 70 ' s. Heckl [2] consi- 
dered a system of beams coupled together to form a 
regular grillage and demonstrated that waves propagate 
only in some frequency bands. Mead [3] included the 
effects of damping in the wave propagation theory for 
periodic beams and later, discussed the nature of the 
propagating waves and their possible interaction with 
acoustic waves [4]. 

In the literature referred to above, exact 
haimionic solutions have been found for the equations 
of motion of the periodic system. Mead [5] and Abra- 
hamson [e] have discussed an extension of Rayleigh's 
method, and the Raleigh-Ritz technique, to peimiit a 
study of the propagating waves in non-uniform periodic 
structures* The significance of this is important, for 
it reveals the possibility of applying approximate 
methods, such as finite element _ technique, to wave 
propagation in periodic systems, Malcing use of finite 
element method 6 rris and Petyt [7] have considered two 
types of periodic construction, one a periodically 
supported infinite beam and the other a skin-rib 
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stiTOctrare^ to find the phase constants and then natural 
frequencies* ■ : _ : .r :r ^ w ■ ■ 

In the cases of the periodic structures made 
Of anisotropic materials^ it is quite difficult to get 
exact harmonic solution by analytical methods. In all 
such cases, finite element method can advantageously 
be used for finding the propagation constants and the 
natural frequencies.' Malliclc and Mead [s] analysed a 
periodically supported single layer ring analytically 
using wave propagation, . Reddy [9] has analysed dynamic 
behaviour of layered rings on periodic supports by 
analytical ' approach. No material anisotiopy has been 
considered in the above works. 

In the present work, finite elsnent approach 
has been applied to obtain natural frequencies and mode 
shapes of periodically supported beams and layered rings 
making use of wave propagation theory. The finite ele- 
ment program is developed to coirpare the natural fre- 
quencies of periodically supported beams obtained by 
Orris and Petyt [7], Also it has been used to obtain 
the natural frequencies arid mode shapes of two layered 
periodically suppoirted ring with three supports and 
these have' been compared with those presented by Reddy 
[9], Finally layered orthotropic periodic ring on 
three supports has been considered and the natural 



frequencies and mode shapes are obtaineu. 

The wrh is broadly divided into three chapters. 
The second chapter deals with the finite element formu- 
lation of periodically supported layered ring to study 
the dynamic behaviour and the third chapter deals with 
the results and discussion. The results are presented 
for two layered radially supported ring. Both layers 
are considered elastic but of different materials. 

These are compared with the earlier results by Reddy 
[9]. Finally frequencies and mode shapes obtained by 
finite element program are presented for a layered 
orthotropic periodically supported ring. 



CHAPTER -- II 


FINITE ELEMENT FORMULATIO N 
FOR PERIODIC RING STRUCTURE 

2.1 Intro duction s 

Finite element method is a versatile numerical 
technique to analyze complicated structural problems. 

In recent years, it finds use in almost every engineer- 
ing application. Analysis of vibration behaviour of 
structures is one such area where the technique has been 
extensively used to get accurate frequencies and mode 
shapes. 

This chapter "briefly outlines the theory of 
finite element foimaulation of periodically supported 
rings to find the propagation constants, natural fre- 
quencies and the associated iTiode shapes. Displacenent 
formulation of the S-node isoparametric elament used in 
the present work is described. 

Derivation of the eigen value problem expression 
for finding natural freq^aencies of a ring using finite 
element method and the method of solving the eigen value 
expression are described briefly. Finally, application 
of the method to the two 3-ayered ring with isotropic 
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material and three layered ring ^Arith orthotropic mate- 
rial, for finding frequencies and mode shapes is des- 
cribed. 

2.2 Finite Element Formulation to Find the Propagation 
Constants and Fr^equencies of the Ring with Periodic 
Suppo rts s 

Finite element formulation applied to dynamic 
problems is described in several boohs [ 1 O- 12 ]. In 
particular^ the approach to find the propagation cons- 
tants and frequencies of a periodic structure is des- 
cribed by Orris and ?etyt [?]. The finite element for- 
mulation for ring type periodic structure follows a 
similar procedure which is described briefly in the 
following sections. 

2.1.1 Finite Element ""ormulction by Pisplacement 
Method; 

For the finite element displacement method of 
analysis gf the periodically supported layered rings 
(Figs. 2.1 and 2.3).- one bay of the ring structure is 
divided into an arbitrary number of discrete elements 
(Figs. 2.2 and 2,4). It may be noted that., bay is the 
span of the ring between one support and the next 
support. Within each element a polynomial type of 
expression is assumed for the displacement field. This 
displacment field is 'then related to the displacements 
at the nodal points located on the element boundaries. 
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FIG. 2.2 PERIODIC BAY OF TWO LAYERED RING 
SHOWING FINITE ELEMENT MESH . 
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FIG. 2.3 THREE LAYRED ORTHOTROPIC RING 
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G.2.4 PERIODIC BAY OF A THREE LAYERED RING 
WITH FINITE ELEMENT MESH. 



The development of the element stiffnecis and mass matrices 
is described in the following section, 

2.2.2 Displacement Formulation of S-node Isoparametric 
El ement ; 

The unique description of the displacement 
within each element in terms of nodal values and boundary 
points or internal points of the element is the basic 
step in any displacement finite element formulation and 
can be expressed as 

= [n] {a^®H (2.1) 

where • {uj is displacement vector.. 

[n] is the shape function matrix 

{a^^H is tbe nodal displacement vector of 
the element. 

If the shape functions chosen to describe the 
element geom.etry are identical to those used to prescribe 
function (displacement) variation, then the element is 
te3rmed iso-parametric. A typical isoparametric element 
is shown in Fig. (2.5). The functional property of the 
inteirpolation function is that its value in the 

natural coordinate system is unity at node 'i' and zero 
at all other nodes. The formulation given below follows 
general pattern of derivation suggested by Bathe and 
Wilson [ 12 ] and Zienkiewicz Ils]. For a two dimensional 
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FIG. 2.5 ISOPARAMETRIC ELEMENT TRANSFORMATION 
FROM LOCAL TO GLOBAL COORDINATES. 
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S-node element,, a qnadiatic variation of the displace- 
ment and geometry are assumed and interpolation func- 
tions to Hg are given in Appandd-X A„ Thus, it has 
16 displacement degrees of freedom- t/TO degrees of free' 
dom being for each norleo Polar coordinates r, © and 
natural coordinates 3 ,. , t^. era roletad as follows* 


x-i 

li 

u 

^1 

■"a 

^2 

•f 

3 

^3 ■■■ ""4 ^4 '^5 

"^5 



+ 




^3 

^3 


(2.2a) 

e = 

©. + 


©^ 

4* 


-!- N , 0 , + 

©r- 

+ ©^ 
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2 
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3 

3 4 4 5 
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+ 

•^7 
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(2,2b) 


Similarly, the displacements expressed in nodal values 
ares 


u = u^ + U 2 Ng Ug -r- u, u^ ■+ 


.j. Ng ug 


^6 ^6 
(2.3a) 


V = v^ t v^ t Ng V 3 V K, v^, Ng Vg + Ng Vg 




(2.3b) 

Thus, the above two cxnaatSons may be written 
in matrix notation as 

(2,4) 
(2.4a) 


i u^®^J = [n] £ a 


where 


[n] = 


( e) 

i f 


V v*' 


0 O N,, 0 N, O N_ O 0 N„ 0 N_ O 
1 2 3 & 6 7 8 

0 0 0 O 0 Ng O Ng 0 0 Ng 


(2.4b) 
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T 


’^2 ^2 ""3 ^3 ^^4 ^4 ^5 "^5 '-"e ^6 

u^ v^ Ug Vgj (2.4c) 


The element; strain vector for plane stress case is 
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(2. 5b) 


The strains at any point within an element are 
related to the nodal displacements by 


(e) 


= [b] 


(e) 
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£ ® J 

where strain displacement matrix is given by 
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(2.6b) 

For plane stress, the elasticity ( stresses train) 
matrix for isotropic material can be written as 
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FIG. 2.6 LONGITUDINAL AND TRANSVERSE DIRECT- 
IONS OF FIBRES IN LAYER 3 AND 1. 




FIG. 2.7 MATERIAL AXIS L, T, Z FOR A LAMINA 
(ORTHOTROPIC LAYER) 



18 
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E 
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■ V 
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V 

1 

0 


O ! 

! 
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(2.7) 


wnere 


'B' is the Young's modulus of elasticity 
V is the Poisson's ratio 


For ortho tropic material, the stress-strain 
matrix for plane stress case when the fibres are oriented 
as shov/n in Fig, (2.6) can be x^rritten as 


where 




^11 ^12 ° 
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i 

d/ 

1’ 
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^21 ^22 ° 

1 (2.8) 
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TL 


Q 

Q 

Q 

Q 

Q 

is the Young' s rrodulus in the direction of fibre 
axis ' 

is the Your^'s modulus in the perpendicular 
direction of the fibre axis 
is the shear modulus 
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is the Poisson ratio in the direction of the 
fibre 

is the Poisson ratio in the direction transverse 
to the fibre axis 


For orthotropic material, the stress-strain 
matrix for plane stress case when the fibres are oriented 
as shovm in Fig, (2,7) can be written as 


[d,] = 


1— V 


ZT 


ET 


O 


V 

1 

O 


ZT 


O 

0 

1-V. 


(2,9) 


ZT ,! 


where is the Young's modulus in the perpendi- 

cular direction of fibre axis 
and Vgrp is the Poisson ratio in the Z direction 


Fig, (2.7a) 

(e) 

The stiffness matrix K and mass matrix M 
element can be evaluated as 

s: / b'^ D B cJ/ 


(e) 


M 


Ce) 


I P n’^ N dV 


of the 

( 2 . 10 ) 

( 2 . 11 ) 


where is density of the material of the elenent. 

The matrices in general are evaluated by 

numerical integration, using Gaussian quadrature formula. 

For plane stress condition, thicTcness 't' being constant, 
(e) (e) 

the matrices K and M can be obtained as 



20 


(e^ 

= t E S a, oc. F. . (2.12a) 

i = l j:^l ^ J 


where 

P. . 
3-3 

T 

= B, . D B. . r. . Det J. . 

3-3 3.J ^ IJ 

(2.12b) 

and 


n n 

:=t S E a. oc.Pf. 

i = l 3. J XJ 

(2.13a) 

where 

^ij 

= P [Njj]’’ Det . 

(2. 13b) 


In the above equations/ is the radial 

J. J 

4“V\ 

coordinate of ip Gaussian point/ and n are the number 
of integration points in each coordinate direction/ 
Det.J^j is the determinant of tlie Jacobian matrix of 
transformation, 

are Gaussian quadrature x\reights of the 
i^ Gaussian point, 

F^j / P^j are function values to be evaluated at 
Gaussian points z 

For evaluation of the stiffness matrix of the 
element, a 2 x 2 numerical integration is used. For the 
evaluation of the mass matrix of the element, a 4 x 4 
integration order leading to the summation of functions 
evaluated at 16 sait^ling points is used, 

2,2.3 Development of Assembled Stiffness and Mass 
Matrices of the Structures 

The assembly of the stiffness and mass matrices 
of all the elements is carried out to get the cpmplete 
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STIFNESS MATRIX 





MASS MATRIX 


BEFORE applying BOUNDARY CONDITIONS 



STIFFNESS MATRIX MASS MATRIX 

AFTER applying BOUNDARY CONDITIONS 


FIG. 2.8 APPLYING BOUNDARY CONDITIONS FOR 
STIFNESS AND MASS MATRICES. 
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• hji, , 

•''^titfness and mass matrices of one bay of the ring. If 
t^ector of nodal displacorents for the complete ring 

?•'* rt 5* ^ 1 L.1 1 ' ■ /- 

s a J then the nodal displacements vector of an 

ndividual element { is related by a transforma- 

*'™^icniori. matriec [t_. ] given by 


f V e>^ 

i a-, s 




(2.14) 


^^'^x^^7'h er* €3 


(2.14a) 


•«»w«i*-wher*e 'n' indicates total number of nodes for one bay 
aBao**of tfie ring. 


The stiffness matrix, for one complete bay of 
the ariiig' consisting of "NE^ elements is then given by 


NE / ^ 

[k=] = E [T.P [k‘®>] [tJ 
i=l 


(2.15) 


A similar expression gives the mass matrix 


j| [m*^] of one conplete bay of the ring. 
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2-2.4- Boundary Conditions ° 

After tlie assembly of stiffness and rrvass 
matrices# displacssnent boundary conditions are imposed 
on them, thus making them non-singular. In this section, 
the application of different boundary conditions for 
one bay of a two layered and a three layered ring are 
described. There are different ways of imposing the 
boundary conditions. The method used in the present 
work is described as follows. 

"tin 

Consider i^ degree of freedom to be constrained. 
The procedure is to make all elements of the stiffness 
and mass matrices corresponding to i row and column as 
zeros except the diagonal element of i"^^ row. This dia- 
gonal element in the present case is chosen to be suffi- 
ciently high and corresponding element in the mass matrix 

is made unity. This ensures, a unit eigen value corres- 

til 

ponding to i degree of freedom to be eliminated in 
finding fundamental frequencies by subspace Iteration 
technique of solving eigen value problem. 

The boundary conditions for monocoupled system 
are applied to the node on the support boundary at mid- 
plane radius. Both the radiad. and tangential degrees 
of freedom are constrained. 
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For Bicoupled case, only radial degree of 
freedom is constrained at the node on the bDundary 
lying on midplane radius of base layer 2, 

For a three layered monocouplad system, the 
boundary conditions are applied in a similar manner, 
only difference being that, the node on the boundaries 
on the midplane radius of base layer 3 are considered 
instead of base layer 2 as in the case of monocoupled 
case. 

2.2,5 Formulation of Eigen Value Expression for 
Propagation Constants and Eigen Values" 

Considering the periodic ring structure shown 
schematically Fig, (2.9), each periodic component AA'BB'', 
could be described by an identical arrange- 
ment of elements. These periodically occurring element 
grids are connected by the displaconents on the boun- 
daries AA' , BB', , 

The assembled vector of nodal displacanents 
for one of the sections, AA'BB' will contain degrees of 
freedom corresponding tp. nodes on the left' hand boundary 
{ a^i , the right hand boundary £ aj^V , , and all other 
nodes, { a^j , Thus, the vector of nodal displacements 
for the section after the application of constraints is 
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31 27 


NUMBER OF NODAL 
DEGREES OF FREEDOM 
IN INTERMEDIATE 
PORTION 76 




14 10 


UMBER OF NODAL 5 
EGREES OF FREEDOM 
N R BOUNDARY 10 


12 0 ” 


NUMER OF NODAL 
DEGREES OF FREEDOM 
ON L BOUNDARY 10 


FIG.2.9 LOCATION OF NODES ON L AND R 
BOUNDARIES. 
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i i 


= i a-. 


(c) 


(c ) 


(2.16) 


Corresponding to each nodal displacement^ this will be 
a generalised force. The force vector for .the periodic 
section may be defined as 


(2.17) 

Whan the periodic systan vibrates harmonically 
with angular frequency o) , the equation of motion is 

([k^^)] -0)2 = {F J (2.18) 

where and are defined with equation_ (2,15)/ 

c4^^3by equation (2,14a) and t F I by equation (2.17). 

At this stage/ in order to facilitate the 

general solution of equation (2.13)/ the displacenents 
(c ) 

{ a I and forces (F) are allowed to talce conplex values. 

When a free wave propagates through an infinite 
structure/ £F^i is zero. The ratio between correspond- 
ing displacements in adjacent bays of the structure is 
equal to ^ z where is the coirplex propagation cons- 
tant. The nodal forces at the exterior nodes of the 
section, } and {Fj^) are rot zero, since these 


{P 5 = 


F, 
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are internal forces in che over all structure. 

The nodal displacements along the left hand 
boundary and those at the corresponding right hand 
boundairy of the section are related by the propagation 
constant. Thus, if the node numbering sequence on the 
right hand boundary is identical to that on the left 
hand boundary, the nodal displacements are related by, 

la^^H = e^ (2.19) 

Similarly, it may be shown that for equilibrium between 
adjacent sections the nodal forces and moments at the 
left hand and right hand boundaries are related by, 

( 2 . 20 ) 

vjhen equations (2,19) and (2.20) are substituted into 
equation (2,18), with { = 0, the e_uation of motion 

of the periodic section becomes 

([k^c)] > 0)2 

The assenbled mass and stiffness matrices, 

[k*^] and may be partitioned bo correspond with the 

three types of nodes, those on the left hand boundary, 
those on the right hand boundary, and all other, inte- 
rior, nodes. Thus, 


^ i 

i 

r p 

; L 

i 

{ 

1 

o 

n 


■ u 


( 2 . 21 ) 
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[k=] = 


[m=] 



^<LI 

\r 

^IL 

H 

H 

^R 

_^RL 


Kr,r. 

RR 


^\l 

^'tR 

M 

^TL 


M 

IR 

RL 

M 

RI 

M 

RR 

( 2 • 19 ) 

and ( 2 

. 20 ) ■ 


and 


( 2 . 22 ) 


(2.23) 


of degrees of freedom present for the constrained system 
(2.21) may now be reduced by the number of degrees of 
freedom along the right hand boundary of the section. 

The matrix equation then found is 


5 K. 


.1^ 


] 


”lR + »RL j “lI hg ] 


II 


(c) 

( c) i 


(c) 


“IL + “IR 


I a. 


(c) ! 


or, in abbreviated form, 

[K(jL(,)] - w2 [M(ja)] 


i Jc) j 

i ^ 

I _ ( c) 

I^I j 


o 


(2.24) 


(2.2 5) 


where [k(^i), [M(]a)] are real functions of tlie complex 


variable jLt 
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ThG set of equations ( 2 *2 5 ) represent an eigen 
value probleiTu Corresponding to any value of u there 
will be a discrete range of frequencies for which a 
vrave motion X’/xth that particular propagation constant 
is possible-. 

Relation (2,2 5) may be rearranged to create 
an eigen value problem for jii, given the frequency w. 

Mead [ 4 ] has demonstrated how a rearrangement of this 
type is possible, Hovrever for the ring under considera- 
tion, the propagation constants are obtained by a much 
simple procedure given by Mallih and Mead (8). 

For the ring, the special type of harmonic 
waves propagate without any reflection as the structure 
does not have any ends. Hence, the ring can be consi- 
dered as an infinite periodic structuie except that the 
propagating v-uves should have the same phase and magni- 
tude after tx'avelling once around the complete ring. 
Consequently, frequencies (natural) at which harmonic 
waves can propagate are not given by the entire propaga- 
tion bard but are confined to a discrete set given by 

= 0 , = + 2j7r/N = 0 , 1 , 2 ,..,N) 

.... (2.26) 

where N is the nxamber of bays; in the ring and the sub- 
script n refers to the value at a natural frequency. 
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As seen from rhe above expression (2,26)^ the 
propagation constants always occur in pairs (positive and 
negat.we) and the number of pairs are equal to the 
minimum, number of coupling coordinates between the adja- 
cent periodic elements (ref. 5). 

Orris and Petyt [?] have pointed out hov? tedious 
it would be/ to solve a complex eigen value routine when 
[kCh)] and [M(ia)] are of high order and in particular 
vxhen [mC/i)] is banded rather than diagonal. Since the 
natu,ral frequencies of the finite structure can be deter- 
mined from the variation of the purely imaginary part 
of the propagation constant, it can be noted that, in 
assessing the modal response of the structure, the fre- 
quency variation of the imaginary propagation constant 
is the more important ruantity. If ii is purely imaginary, 
the matrices in equation (2.2 4) and (2.2 5) are Hermitian. 
The solution of the complex eigen value problem can be 
reduced to thcit of a real symmetric matrix, as described 
by Willcinson [is]. 

Let [k(A)] = + IK^], [m(ii)] = [m^ + iM^3 

(2.27) 

/ ^^i^^ = Cb^ + ib^i f 

Then equation 2.2 5 becomes 

([k^ t iK^] + iM^])i ° 


(2.23) 
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By separating real and imaginary parts, and 
combining the two sets of equations, then obtained, the 
matrix expression found is 

o r ! b^ 

- W 2 [ 1 . 

M 

As the matrices [k( 4 )], [m( 4 )] are Heirmitian, 
and [m'^]= Thus, equation (14) 

represents a real symmetric eigen value problem. The 
equation has been replaced by a real matrix expression, 
but at the expense of doubling the size of the matrices 
involved. The latter disadvantage is however offset 
by the fact that, in general more efficient computer 
routines are available to solve real eigen value problems. 

Corresponding to each value of ji in the real 
matrix expression (2.29), there will be a discrete set 
of frequencies that occur in equal pairs. Each frequency 
will have associated with an eigen vector. This eigen 
vector defines the wave motion in periodic section at 
that frequency. 


=0 (2.29) 



2.2.7 Subspace Iteration Method for Solving 
Eigen Value Problem: 

This method is effective for finding the first 
few eigen values and corresponding eigen vectors of 
large eigen value problems whose stiffness and mass 
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matrices have large bane widths. The various steps of 
this method are described briefly in Appendix B, A 
detailed description of the method can be found in Bathe 
and Wilson [lo] and several other books [ll, 14, IB]. 

2.3 Application of. the Method to Two Layered 
Elastic^ Ring With Periodic Supports •• 

A two layered ring, with- different materials 
as shown in Pig, (2,1) is considered. The span between 
supports 1 and 2 is shown separately in Fig. (2,2). This 
is divided into 12 elements. Each of these elements is 
a 8~node isoparametric plane stress element in polar 
coordinates (r,9). Radial supports are provided at the 
midplane radial nodes at the left and right hand boun- 
daries as shown in Fig, (2,2). 

Lifferent support conditions o.re considered 
to get the ocrresponding frequencies and eigen modes, and 
are compared with those obtained by Reddy [ 9 ], 

In the first case both the radial and tangential 
displacements (u and v respectively), of the supports 
are constrained. This results in one coupling coordi- 
nate between the two adjacent bays. Effect of 

increasing the rotational stiffness of tie ring by in- 
cluding two rotational springs at the two boundaries 


is also considered. 



33 


In the second case, only the radial displace- 
ment u is constrained at both 'the supports. This 
results in two coupling coordinates v, ^ between the 
adjacent bays. Since there ex5.sts a rigid body rotation 
in this case, the method described by Zienkiwicz [is] 
is used to overcome the problem of singular stiffness 
matrix. This method primarily involves shifting the 
frequencies by an amount cc, and subtracting the same 
after solving the eigen value expression. 

For all the cases considered above, mode 
shapes have been plotted. 


2.4 Application of the Method to Three layered 
Ortho tropi.c Ring with Periodic Supports s 

A three layered ring as shown in Fig, (2,3) 
has been considered for the present wor'-:, with outer 
and inner layers having the fibre orientation as shown 
in Fig, (2.6) and the middle layer having the fibre 
orientation as shown in Fig. (2.7), 

One bay of the three layered ring with three 
supports is considered for obtaining the natural fre- 
quencies and mode shapes using ^wave approach'. Each 
one of the layers, in the 3 layered bay, is divided into 
six '8— node isoparametric elements, to>-alling upto 
13 elements. Since each of the outer and inner layers 
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has different fibre orientation, it is interesting to 
see how the 'three layered ring behaves when a free wave 
propagates. Frequencies and mode shapes are obtained for 
monocoupling case where both radial and tangential 
displacements of the nodes situated on -the midplane 
radius of Iciyer three are constrained giving rise to 
only one coupling coordinate §§ . 



CHAPTER - II 


REjSUL T S AND PI SOJ SS ION 

3.1 Introduction 

The accuracy with which the propagation cons- 
tants and their corresponding natural freguencies can 
be obtained by finite element method is first cheched 
by comparing the results for periodically supported 
beam against tlie results of [?] and other \>rorkers«Next the 
reisults are conipared for two cases of a periodically 
sxipported two layered elastic ring with those obtained 
by [ 9 ]. Finally# frequencies and mode shapes of a three 
layered orthotropic ring obtained by the present method 

are sliown, , 

i 

3.2 Preguencios for Periodically Supported Beams 
Rof, to a Pig* ( 3. 50 )of Beam: 

To check the accuracy of the developed PEI4 
computer program# the purely imaginary part of the 
propagation constant is plotted as a function of non- 
dimensional frdxjuency. The curve Figure C3»5y is 

^^?ork, 

found to be exactly matching with that presented in earlier^ 
Pctyt's results need correction by a factor of 40. 
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?.? Two-Layered Elastic Ring on Three Periodic 

s-upports ” 

The present, finite element analysis is used 
to obtain tlie natural frequencies and mode shapes of 
tvjo layered elastic ring on three periodic radial 
supports as shovm in Pig. (2.1). The numerical data 
presented by Reddy [ 9 } is "coh-sidered hare for better 
comparison of results.. This' is presented as follox-^s: 

N = 3 (three supports), = 0 and 2* 
whQ.re ratio of half thickness to 

mldplanc radius of layer 2 , E^ and E^ are Young's 
modulii of elasticity for outer and inner layers, and 
p 2 ate mass densities corresponding to layer 1 and 
layer 2 , N.=s 3 (number of supports of the ring)^ 
is a nondiraensional rotational stiffness applied in 
tho form of spring at both the boundaries of one bay 
of the ring* . . 

3.3,1 Mo no coupled Two Layer Ring on Tliree Supports 

with 5^ = 0 

. . ''i' 

The nondimensional natural 'frequencies obtained 

are tabulated in Table 3.1 along wi.th those obtained by 
Reddy [ 9 ] for the phase constants = 0, 27 t/ 3,’ where 
the nondimensional natural frequency is given by 

^ a •/ . 
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where 


m, 


2 is mass oi -the base layer 2 over a length 


where 


of R. . 

1 

= 2 t hj P 2 Rj 

6 2 is density of material of base layer 2 
b is width of the ring 
^2 is the bending stiffness 


2 

3 


Ej b 


where ^ 




h 2 is bay thickness of base layer 2 

E 2 is Young's modulus elasticity of base 
layer 2 

R 2 is radius of the ring. 


Prom the Table 3,1 it can be observed that the 
error by the present analysis is within 0 , 5 ? of those 
presented by Ref. [ 9 ]. Also for the positive going 
wave (ji.„ = -27 t/ 3) and negative going wave (%-, = +2Tt/3) 
the frequencies are found to be exactly matching. The 
number of subspace iteration vectors used to obtain 
first four frequencies accurately are 8 and this has 
taken 16 iterations to converge. The computer process- 
ing time is noted to be 480 sec, (computer time) mode 
shapess \ 

The eigen value expression (Eq. 2,29) derived 
in sec. 2,2,5/ gives complex eigenvectors £ i 
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for the same freqoaency tor a given value of phase cons- 
tant. This is because of the doubling the size of the 

eigen value expression given by Sg. 2.2 5. We consider 
b-^ 

only ^ i as this is the vector reguired by expression 

2.2 5- For = + 27 t/ 3 and n. = - 2n/3 also, vectors 

i i+i i -iJ are obtained (Table 3.5). 

b b^ 

Tt is observed from the results obtained that, 

in the present analysis, the eigen value , at ju, = +27 t/3 

in 

and ^ = -2Tr/3, is not repeated. But it has been estab- 
lished Ref. [s] that there are degeneracies at these 
frequencies. 

In the •work by Reddy [ 9 ], it is noted ■t±iat for 
a standing mode to exist at a natural frequency the 
moments M and M of the positive goir^ wave and nega- 
tive going at both the supports should be either in 
phase or anti phase with each other. This is earlier 
been verified experimentally also by Mallik and Mead [s]. 
Since there are tw conditions at which the mode can 
exist at 0, stated as above, there is bound to be 

a degeneracy of modes at the frequencies at which 

The present finite element analysis does not 
directly lead us to the second eigen vector at these 
degenerate frequencies. Rather knowing that a second 
eigen vector at the same frequency for 0 may exist. 
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and also that the standing waves are sought for/ an 
explanation for tne existence of second eigen vector 
is given in the following. 

The total deflection vector vjill be a sum of 
positive and negative going waves for the first bay of 
the ring, given byf bj = {{ + i£ b^fj +{ { b^j 

+ i£ b^ / whereCb^"^ is real part of = -27r/3 and 

i b^"*) is the imaginary part of the positive going wave 

ZT*" zr 

andjb ^ anc^ b ^ are the real and imaginary* parts of 
negative going wave at = +27 t/ 3. For standing waves 
to exist, the phase difference between positive going 
wave and negative going x-jave has to be aero leading to 
either ^b^'j = -ib^“i ori b^"^ = -{IcF}. 

Since the results presented in Table 3.5 show 
that { + ii b^'^j is not the complex conjugate of 

eithori { b’^'i + or ~ { Cb^l + i{b^“)j . It is 

guessed that at . = 2tt/ 3, there is another eigenvector 
{ Ib^*^!- iCb^’^B and similarly at = - 2n/3 there is 
another eigen vector equal to { £b“ i- iCb i} . With 
this expectation, the standing wave at h = - 27 t/ 3 is 
given by ib^'^i and at y. = + 27 t/ 3 by C b"^> . The 
results are plotted in Fig. (3.1) for each bay of the 
ring and these are found to be matching with shapes 
presented by Reddy [ 9 ]. This supports the eJ^lanation 


given above 
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TABLE rt.l 

Natural froquencies of two layered elastic ring 
with mono coupling and = 0 

= 0*2 , E^/E^ = 0.001, 6 ^ =: 0.02 

\/^2 “2/ N =: 3 (Number of supports) 


Mode 

I 

^-in 

5 Q 

5 n 


No. 

§ 

5 

5 

t 

1 n 

! Finite 1 

jj el ernent 5 

X results 5 

! ^ 

Ref. (9) 
results 

1 


0 

5.9 5? 

5.973 

2 


+27T/3 

3.947 

8.976 

3 


+27T/3 

16.412 

16.331 

4 


0 

20.210 

20,13 
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TABLE 3.2 

Natural frfxjuencies of two layered elastic ring 
with monocoi-pling and = 2 

6/62 = 2, N = 3 (Number of supports) 



1 

0 

6. 602 

6.619 

2 

+27T/3 

9 . 148 

9 . 179 

3 

+27T/3 

16,788 

16.750 


0 20.21 20.13 





42 


TABLE 

Natural frequencies of two layered elastic riiig 
with Bi coupling and = 0 

P/P2 = = 0*^01/ ^2 = 

^2 “ 2 / N = 3 (Number of supports) 


Mode 5 % 

No. J ^ ^'Finite T Ref .(9) 

jj 5 element i Results 

5 I results ^ 


1 0 

2 + 27T/3 

3 0 
+ 27t/3 


1.416 1.414 

6,570 6.581 

9.002 9.0 31 


4 
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TABLE :^.4. 


Natural frequencies of three layered orthotropic 
ring with monocoupling and = 0 (same material for 
all 3 layers) 

^3 = 0.015 , ~ = ^ / N = 3 (Number of 

supports) 


Mode 

No. 

ILL. 

in 

Q 

n 

Finite element 
results 

1 

0 

23.936 

2 

+ 27T/3 

37^668 

3 

+ 2 Tt/ 3 

59.971 

4 

0 

74.031 
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TABLE e; 

DISPLACEMENT VECTORS FOR THE TV70 LAYERED RING VTITH 

MONO COUPLING AND K = 0 

r 

Case (i) = 0 Frequency? 5.9.53 


Node No. First Bay Second Bay Third Bay 


2 

0.0 

0.0 

0. 

.0 

10 

-0.4 

-0.4 

-0, 

.4 

18 

-0.4 

-0.4 

-0, 

.4 

26 

0.0 

0.0 

0. 

.0 

34 

0.4 

0.4 

0. 

.4 

42 

0.4 

0.4 

0, 

.4 

50 

0.0 

. 0.0 

0. 

.0 


Case (ii) a^^=27T/3 


-*XUv\ 

Frequency ? 8.9 47 


2 

0.0 

0.0 

O.O 

io 

-0.09 

-0.61 

0.71 

13 

-0,45 

-0.48 

0.93 

26 

-0.28 

0.24 

0,03 

34 

0.26 

0.64 

-0.91 

42 

0.47 

0.28 

-0.76 


0,0 

0 . 0 < 

0.0 


Case (iii) 

‘‘ln= +2n/3 

Frequency : 8,947 

2 

0.0 

0.0 

0,0 

10 

-0.49 

-0.26 

0.75 

18 

-0.32 

40.02 

0,80 

26 

-0.17 

40.36 

-0.12 

34 

0.71 

40.17 

-0.83 

42 

0.72 

-0.14 

-0. 58 

50 

0.0 

0.0 

0.0 
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TABljS 3.6 


DISPLACEMENT VECTORS FOR THE TWO LAYERED RING VHTH 

MONO COUPLING AND K = o 

r 

Case Li) = 0^ Frequency = 6.602 


Node No. 

First Bay 

Second 

Bay Third Bay 

■ 2 

0.0 

0.0 

0.0 

10 

0.84 

0.84 

0.84 

18 

0.87 

0.87 

0.37 

26 

0.0'" 

0.0* 

0.0" 

34 

-0.87 

-0.87 

-0.37 

42 

-0.34 

-0,84 

■ -0.84 

SO 

0.0 

0.0 

0.0 

Case (ii) ti. = 

+ 271/3 

%rLv\ 

Frequency = 9,148 


2 

0.0 

0.0 

0.0 

10 

- 0 . 50 

0.74 

-0.24 

18 

-0.38 

0.95 

-0.57 

26 

+0.23 

O.Oi- 

-0,22 

34 

■10.52 

-0.9 5 

-^.43 

42 

+0.20 

-0.73 

0.54 

50 

0.0 

0.0 

0.0 

^Case (iiii* M-. 

= - 2.71/3 

wTLv^ 

Frequency = 9,148 


2 

0.0 

0.0 

0.0 

10 

-0.07 

0,69 

-0.62 

18 

0,22 

0.69 

-0.90 

26 

0,26 

-0,16 

-0.09 

34 

-O.O'S 

-0.80 

0.85 

42 

-0.27 

-0.50 

0,7 5 

50 

0.0 

0.0 

0 .0 
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TABLE ?,.l 


DI.SPLACMENT VECTORS FOR THE TWO LAYERED RING WITH 


BI -COUPLING AND K •’ = 0 

r 

Case (i) II. = 0 
in 


F requency = 6 . 57 


0) 

o 

No. First Bay 

Second Bay Third Bay 

2 

0.0 

0.0 

0.0 

10 

-0.92 

-0.92 

-0.92 

18 

-0.92 

-0.92 

-0.92 

26 

0.0 

0.0 

0.0 

34 

0.92 

0.92 

0.92 

42 

0.92 

0.92 

0.92 

50 

0.0 

0.0 

0,0 


Case (ii) = 




+ 271/3 

Frequency = 1.41 

2 

.O.O 

0.0 

0.0 

10 

-0.26 

-0.16 

0.43 

18 

-0.35 

-0.4i 

0.76 

26 

-0.29 

-0.57 

0.37 

34 

-0.14 

-0. 57 

0.72 

42 

-0,12 

-0.37 

0.38 

50 

0.0 

0.0 

0.0 


Case (iii) U- 

= -271/3 

— A-v> 

Frequency = 1.41 


2 

0.0 

0 . 0 , 

0 

.0 

10 

■ 0.26 

^ 0^16 

40 

.43 

IS 

0.36 

- 0.30 

4O 

.76 

26 

0 .:^ 

■ - 0 . 57 

r -0 

.002 

34 

0,15 

- 0 . 57 

-0 

.72 

42 

0.01 

- 0 . 36 

-0 

.38 

50 

0.0 

0 .0 

0 

.0 
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TABLE 

DI SPLACEMEI-TT VECTORS FOR THE THREE LAYERED ORTtlOTROPIC 

RING V/IIH iMONO COUPLING AND K = 0 

r 

Case (i) 0 PrequencY^ 2 3.98 


®ode 

No. 

First Bay 

Second Bay 

Third Bay 

' 2 


0.0' 

0,0' 

0.0' 

13 


0.37 

0.37 

0.37 

24 


0.36 

0,36 

0.36 

35 


0.0 

0.0 

0.0 

46 


-0.36 

-0.36 

-0.36 

57 


-0.37 

-0.37 

-0.37 

68 


0.0 

0,0 

0.0 


Case 

( ii ) = 

. — iX-v\ 

+2n/3 Frequency = 37,66 

2 


0.0' 

0.0 

0.0 

13 


0.31 

0.34 

0.03 

24 


0.22 

0.47 

0.2 5 

35 


-0.18 

0.16 

0.23 

46 


-0.40 

0,35 

-0.02 

57 


-0.18 

0, 38 

-0.19 

68 


0.0 

0.0 

0.0 


Case 

( iii ) = 

- 27T/3 Frequences 37.66 

2 


0.0 

0.0 

0.0 

13 


0. 36 

-0.08 

-0.25 

24 


0.46 

-0,2. 

-0.15 

35 


0.15 

-0.25 

0,20 

46 


•^0.4 

40.1 

0.35 

57 


-0.4 

0.24 

0.13 

68 


0.0 

0.0 

0.0 

















.3.5(b) VARIATION OF IMAGINARY PART OF PROPAGATION 
rnwqTANT AND NONDIMENSIONAL FREQUENCY. 
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3.3.2 Mono coupled Two Layer Ring on Three 

Supports with K = 2 ° 

r 

The rotation stiffness is included in this case 
by means of springs at the two boundaries. The natural 
f reguGiicies 3-cr this case are tabulated in Table 3.2. 

And uho Table 3.6 shows the displacoment vectors for the 
three bays of the ring . The mode shapes are plotted as 
shown in Fig, 3.2 and these are found to be matching well 
with those obtained by Reddy [ 9 ]'. 

The frequencies in this case also are found to 
be with 0 . 5 ^ accuracy compared to (Ref, 9). The initial 
eigen vector taken for first 4 frequencies to converge 
are eight, and the CPU time in this case also is found 
to be 480 sec, for a total of 15 subspace iterations, 

3.3.3 Bi -Coupled Ta-o Layer Ring on 3 Supports with K^=0 

The toundary conditions in this case are to 
constrain radial degree of freedom at nodes 2 and 50 as 
shown in Pig, 2.2, Here there is a rigid body rotation. 
So the first frequency obtained should be zero. The 
rigid body mode i.s suppressed as escplained in Section 

2.3 

Tho natural frequencies are tabulated in Table 

3.3 and the Table 3.7 shows the displacement vectors for 
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the three bays of the ring. Prom Table 3.3 it is observed 
that due to numerical error^ the frequency for rigid 
body mode is snown as All the frequencies are found 

to be within 0 . 5 ^ accuracy compared to those in (Ref .9). 
The number of initial eigen vectors chosen in this case 
are 8 to obtain first 4 frequencies ^accurately and this 
has taken 439 sec. of CPU time for 16 iterations of 
convergence. The mode shapes plotted are as shoxim in 
Fig, 3,3 and these also match well with those obtained 
by Reddy [ 9 ], 

3, 5 Conclusions 

Many researchers have found exact harmonic 
solutions for finding propagation constants and fre- 
quencies in th.G case of periodically supported struc- 
tures/ when the material is isotropic. It becomes 
difficult to got harmonic solutions for complex problems 
analytically when the material is anisotropic. The 
finite el emont method presented in this work has revealed 
that the material orthotropy can easily be included in 
the analysis to find natural frequencies and correspond- ' 
ing mode vshapes. 
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^P 5NDIX h 

.SHAPE FUNCTIONS OP B-NOBE ISOPARAMETRIi 

“1 = 4 - (|) Hg 

”3 = 4 <1-U’ ■ ^2’ '’5 - ‘2^ **6 

N3 = I (l-.») (1+tj) - ci) Kg - (1) K3 

= i ( 1 -s) (l+t^) - (i) Kg - ( 1 ) Kg 

= 1 (I-S^) 

Ng = 2 ^ l“t|) (l+s) 

= ~ (l“S^) 

^8 “ 2 


ELEI'lENT 
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APPENDIX B 
RAYIjEIGH-RITZ 
SUBSPACS ITSPATION METHOD 

-i-his method oi: solving eigen value express 
of the form 

{ [a] - co^ [b]J I X> = 0 

where [a] and [b] are real symmetric consists of the 
f o 1 lowi ng steps » 

Step 1 : Start with g initial iteration vectors 

^2 ' " ■> ^ P / where p is the 

number of eigen values and eigen vectors to 
be calculated. Bathe and Wilson suggested 
e value of c. = min (2p,p+8' for good con- 
vercTence, Define the initial modal matrix 
as 

X^J = [ X2,..„X^J 

and set the iteration number b = 0. 

Step 2 5 The following subspa.ee it -ration procedure 

is used to generate an improved modal matrix 

(a) Find [ J the relation 

[a] [ w [x^] 



59 


(b) 


(c) 


(d) 


Step 3" 


Compute 

“ I- b] [ and 

[B] [ 

Solve for the eigen values and eigen vectors 
of the reduced system 

= [\+i3 

and obtain and 

Find an improved apparoximation to the eigen 
vectors of the original system as 




If and denote the approximations 

to the i^^ eigen value in the iteration 3c-l and 
1< respectively^ we assume convergence of the 
process whenever the following criteria are 


sati sf ieds 


v(h+l) _ 

'^i ” ^^i 


, (k+1) 

h 


^ Q'pSf 1 = 1/2^ ♦• •p 


where eps = 10“^. It is to be noted that 
although the iteratibn is performed with 
q vectors (q >P )/ the convergence is measured 
only on the approximations predicted for the p 
smallest eigen values. 
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